Stripes in the Ising Limit 
of Models for the Cuprates 
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The hole-doped standard and extended t-J models on ladders with anisotropic Heisenberg interac- 
tions are studied computationally in the interval 0.0 < A < 1.0 (A = 0, Ising; A = 1, Heisenberg). 
It is shown that the approximately half-doped stripes recently discussed at A = 1 survive in the 
anisotropic case (A<1.0), particularly in the "extended" model. Due to the absence of spin fluctua- 
tions in the Ising limit and working in the rung basis, a simple picture emerges in which the stripe 
structure can be mostly constructed from the solution of the t-J model on chains. A comparison of 
results in the range 0.0 < A < 1.0 suggests that this picture is valid up to the Heisenberg limit. 

PACS numbers: 74.20.-z, 74.20.Mn, 75.25.Dw 
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In recent years, evidence is accumulating that at least 
in one family of high-temperature superconducting com- 
pounds (La2_ x Sr x Cu04) one-dimensional (ID) charged 
stripes are formed upon hole doping of the insulating 
parent material pi] . Although the presence of stripes in 
other compounds such as YBa2Cu30g+5 is still contro- 
versial M, a large theoretical effort has focussed on the 
search for stripes in models for the cuprates. Early re- 
sults reported stripes in Hartree-Fock treatments of the 
Hubbard model and also in the phase-separated regime 
of the t-J model upon the introduction of long-range 
Coulomb interactions Hf|. However, these stripes have 
a hole density nh=1.0 in contradiction with the nh=0.5 
density found experimentally |l|,^|. Improving upon this 
situation, recent studies of the t-J model ||[7| reported 
stripes with nh<1.0 at intermediate values of J/t, some- 
times described as a condensation of d-wave pairs ||. 
However, the mechanism leading to nh<1.0 stripes, and 
even the presence of a striped ground-state in the t-J 
model at intermediate couplings, is still under discussion 
|1. In addition, evidence is accumulating that the pure 
t-J model is not sufficient for the cuprates, and its "ex- 
tended" version with hopping beyond nearest-neighbor 
sites pO|Jll|l is needed to explain PES results for the in- 
sulators 

Very recently, indications of half-doped stripes have 
been found by Martins et al. in the extended t-J model 
p3| . They were also observed at small J/t in the stan- 
dard t-J model, in both cases in regimes where two holes 
do not form bound states. This led to a novel ratio- 
nalization of stripes as the natural way in which spin- 
charge separation is achieved in two-dimensional systems 
p3| , similarly as in the phenomenological "holons in a 
row" picture j|. Moreover, stripes appear to emerge di- 
rectly from the one- hole properties of the insulator [jl3| , 
where holes have strong "across the hole" antiferromag- 
netic (AF) correlations pd| , p^ | in their (frustrated) effort 
to achieve spin-charge separation, similarly as it occurs 
in ID systems fla]. 



In spite of this progress, more work is needed to un- 
derstand these complex striped states. An important 
issue is the role played by fluctuations in the spin sec- 
tor. While the presence of a spin tendency to form an 
AF background is crucial for stripe formation, it is un- 
known whether the fine details of the spin sector (such 
as the presence of low-energy excitations) are important 
for its stabilization. To address this question, here a 
computational study is reported where the spin inter- 
action contains an Ising anisotropy. Our main result is 
that stripes survive the introduction of this anisotropy, 
and a fully SU(2) symmetric interaction is not needed for 
stripe formation. This result is in agreement with recent 
retraceable-path calculations for the t- 3 Z model [[L6| . 

The anisotropic extended t-J model is defined as 



H = J^[S z i S z j + -(S+ i S- 

(ij> 



1 , 

jninj] 



tim(clc m + h.C.), 



(1) 



where ti m is t(=l) for nearest-neighbors (NN) hopping 
between sites i and m, t' for next NN, t" for next-next 
NN, and zero otherwise. The anisotropy in the spin 
sector is controled by A (A = 0, Ising; A = 1, Heisen- 
berg). The rest of the notation is standard. t'=-0.35 and 
t"=0.25 are believed to be relevant to explain PES data 
pp| The computational work is carried out using the 
Density Matrix Renormalization Group (DMRG) [p|Jl7[| 
and Lanczos techniques, as well as an algorithm us- 
ing a small fraction of the ladder rung-basis (optimized 
reduced-basis approximation, or ORB A |l9| ). Results 
are presented in (i) the small J/t region with t'=t"=0.0, 
and (ii) small and intermediate J/t with nonzero t' and 
t", in both cases working at A=0.0, 0.25, 0.50 and 1.00. 
These two regions (i) and (ii) have similar physics Jll| , 
and the extra hoppings are expected to avoid phase sep- 
aration ||,|2(|. 
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To start our investigation, let us analyze the one-hole 
properties of model Eq.(l). Previous studies found a ro- 
bust AF correlation between spins across-the-hole (Cah) 
0,0, namely between the spins separated by two lat- 
tice spacings located on both sides of a hole (working in 
the reference frame of the latter). This curious feature 
was interpreted as a short distance tendency toward spin- 
charge separation JjllJ , and it is believed to be crucial for 
stripe formation ]13| . It is important to clarify if similar 
correlations are still present in the anisotropic case. For 
this purpose, here 4x4 and 4x6 clusters were used with 
periodic boundary conditions (PBC), as well as cilindri- 
cal boundary conditions (CBC), (open boundary condi- 
tions (OBC) along legs and PBC along rungs Q), for sev- 
eral different couplings l^jj] . The results for Cah in the 
one-hole lowest energy state at several momenta are in 
Fig. la using a PBC 4x6 cluster, at a representative cou- 
pling set. At A=1.0 all the correlations are negative, com- 
patible with the results in Refs. jll],[l3). As A is reduced, 
the magnitude of Cah decreases but it remains negative 
for the momenta of most relevance, namely (7r/3,7r/2) 
(the closest on the 4x6 cluster to (w/2,tt/2)) and (0, tt)- 
(ir, 0). In Fig. lb the AF spin correlations are shown at 
(0, 7r) and A=0.0. There is a good agreement with the 
results observed in the A=1.0 limit pM. Similar conclu- 
sions are reached in the (ir, ir) case (Fig.lc) at moder- 
ate anisotropy, and for (tt/3, 7t/2) (not shown) although 
in this case the correlation along the short direction is 
approximately zero. The across-the-hole correlations re- 
main even with Ising anisotropy. 
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FIG. 1. (a) Exact Diagonalization calculations of 
across-the-hole spin-spin correlations (Cah) vs. A for several 
momenta on a 4x6 cluster with 1 hole, PBC, J=0.2, t'=-0.35, 
and t"=0.25. The results shown are for the long direction, (b) 
AF spin-spin correlations for (a) at (7r,0) and A=0.0. (c) AF 
spin-spin correlations for (a) at (n,Tr) and A=0.5. In (b) and 
(c) the hole is projected from the ground-state to the site 
shown, and the dark lines represent AF spin-spin correlations 
with a thickness proportional to its absolute value. 

On the other hand, for momenta (0,0) and (7r, 7t) with 
PBC (and also at kj,=0 and ir, with CBC) a transition 
to a ferromagnetic (FM) correlation occurs at small 



A. This tendency is dangerous for the stripe formation, 
which will not occur in a spin polarized background. Cal- 
culating the NN spin-correlations without hole projection, 
i.e. involving all the links of the cluster for several cou- 
plings and 0.0 < A < 1.0, similar trends were observed, 
including the change of sign for A — > 0.0. This suggests 
that the Cah A-dependence at (0,0) and (ir, tt) for the 
one-hole results, can be ascribed to an overall decrease 
of the AF tendency on the clusters studied, which is re- 
placed by a tendency toward ferromagnetism. A compar- 
ison of data at several couplings and lattice sizes supports 
this interpretation. Fortunately, the results are also in 
agreement with the early analysis by Barnes et al. |2^] 
that reported a reduction in the FM tendency upon the 
increase of the cluster size for the Ising limit of the t-J 
model. Thus, it is believed that the FM state found in 
some sectors of the one-hole problem at finite J will dis- 
appear as the clusters grow p3]. In fact, if a compromise 
between couplings and cluster size is followed to avoid the 
FM region, the qualitative "shape" of the A=1.0 one- hole 
wave-function can be preserved as A— >0.0. In conclusion, 
the FM tendency in some subspaces of the one-hole sec- 
tor is expected not to be detrimental to stripe formation, 
and studies with more holes shown below support this 
view. Nevertheless, care must be taken with the stripe- 
FM competition in these systems. 
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FIG. 2. Spin-spin correlations vs A for the most proba- 
ble ground-state hole configuration (stripe) using 4x4 clus- 
ters with 2 holes, PBC and CBC. (a) Spin-spin correlations 
across-the-hole calculated along the stripe; PBC were used in 
both directions. Squares are for J=0.2 and t'=t"=0.0; circles 
for J=0.2, t'=-0.35, and t"=0.25; and triangles for J=0.4, 
t'=-0.35, and t"=0.25. In all three cases k^ky^O. (b) Same 
as (a) but now the spin-spin correlations across-the-hole are 
calculated across the stripe, (c) Same as (a) but now using 
CBC k H =0 (the stripe runs along the PBC direction), (d) 
Same as (b) but now using CBC k M =0. 

Consider now two holes on 4-leg ladders. Similarly as 
in Ref . |lj| , stripes are formed with the two holes mainly 
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located at two lattice spacings along the rung. Then, cor- 
relations along and across the stripe must be considered 
separately. As in Fig. 2, the introduction of a second hole 
on a 4 x 4 cluster using the extended t- J model stabilizes 
robust Cah both along and across the stripe, in PBC and 
CBC, even for values of A where the one-hole system did 
not have a robust AF Cah- For the standard t-J model 
at J=0.2, t'=t"=0.0, also shown in Fig. 2, the situation is 
less clear and in some cases the correlations become FM 
at A=0.0, but in most situations it remains AF. In the 
extended t-J model the stripe tendency is clearly stronger 
than in its standard version p3|. 
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FIG. 3. Rung hole density (n(r)) vs rung index r using 
DMRG (243 states) on a 4x8 cluster with 4 holes, CBC (OBC 
along the direction shown, with invariance under reflections 
assumed). 

Previous studies at A=1.0 showed the coexistence of 
two or more n? l =0.5 stripes at hole density x=l/8 on 4-leg 
ladders as the number of holes increase in a CBC cluster. 
It is important to show that this stripe phase survives the 
decrease of the spin fluctuations. That this is the case can 
be observed in Fig. 3 where DMRG results for the rung 
density (n(r)) vs the rung label r are presented for a 4x8 
cluster with 4 holes (J=0.3, t'=t"=0.0; CBC). Figure 3 
shows that from A=1.0 down to A=0.25 the two n^=0.5 
stripes are virtually unaltered. On the other hand, at 
A=0.0 the FM tendency already discussed in the 2-holcs 
case of Fig. 2 prevent the formation of stripes in the stan- 
dard t-J model. This problem does not occur in the ex- 
tended version with t'<0 and t">0. 

Let us now analyze the spin structure around the stripe 
once the spin fluctuations are fully turned off (A=0.0). In 
Fig. 4a, a stripe configuration with a clear 7r-shift across- 
the-stripe is shown for a CBC 4x6 cluster with 2 holes 
at A=0.0. This dominant hole configuration was pro- 
jected out of the ground state of an ORBA calculation 
where w 10 6 states were kept in the rung basis. The cou- 
plings (J=0.2, t'=0.0, t"=0.25) differ slightly from those 
in Fig. 2, but the results are representative and they are 
similar to those found at A=1.0 [ p5[ . 
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FIG. 4. (a) Spin-spin correlations for 2 mobile holes pro- 
jected at their most probable configuration (circles) in the 2 
holes ORBA ground-state of a 4x6 cluster using CBC (k y =0), 
J=0.2, t'=0.0, t"=0.25, and A=0.0. Full lines indicate AF cor- 
relations (thickness proportional to absolute value); dashed 
lines indicate ferromagnetic correlations, (b) AF spin-spin 
correlations of the CBC 4x4 cluster with 2 holes projected at 
their most probable distance. This result was obtained keep- 
ing just the 8 rung-basis states with the highest weight in the 
ground state (out of a full Hilbert space with 102,960 states). 
Same couplings and momentum as in (a). This approxima- 
tion emphasizes the ID character of the stripe; note also that 
the stripe formed by this procedure is completely decoupled 
from the rest of the cluster. 

It is expected that the suppression of spin fluctuations 
would simplify the description of a ground state like the 
one in Fig. 4a, either in the rung- or Sz-basis, since less 
states are needed to represent the spin background. In- 
deed the combination of CBC with the use of a rung basis 
(along the PBC direction) leads to a fairly simple de- 
scription of such a state. This can be observed in Fig. 4b, 
where the spin-correlations were calculated on a CBC 
4x4 cluster with 2 holes (same coupling as in Fig. 4a) us- 
ing only the eight highest weight rung-basis states out 
of the full ground state of the system, that has a total 
of 102,960 states. The results are very similar to those 
found with the full ground-state and they capture the 
basic physics contained in the full calculation. The pic- 
ture that emerges is the following: the spin-correlations 
along the stripe are maximized, i.e., the two spins on it 
are locked in a singlet, implying that their correlations 
with the rest of the spins vanish. This means that in 
the "snapshot" of the ground-state (Fig. 4b) the stripe 
is disconnected from the rest of the cluster, emphasiz- 
ing its ID character. In fact the only rung-state that 
contains holes in the eight most dominant states kept in 
Fig. 4b is the ground-state of the 2-hole sector of a sim- 
ple 4-site ring calculation. This establishes a strikingly 
simple connection of the stripe problem with a truly ID 
calculation. The "recipe" to construct a good represen- 
tation of the stripe state is to consider the solution of 
half-doped ID chains as a stripe, and antiferromagneti- 
cally couple the rest of the spins of the plane simply as 
if those stripes would be absent (thus generating the n- 
shift across the stripes). This is a natural generalization 
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to two-dimensions (2D) of the ID spin-charge separation 
concept, as emphasized in Ref. [fl3f , where the spin por- 
tion of the wave function is constructed simply ignoring 
the holes In 2D now it is the stripes that are ig- 

nored by the spins not belonging to them in their wave 
function. It is also important to notice that this picture 
appears to hold, in its main aspects, even at A=1.0. 

For completeness, in Fig. 5a ORBA results (with « 
2xl0 6 states kept in the rung-basis) are presented for 
(n(r)) vs r, showing two n^=0.5 stripes with two holes 
each. The calculation was performed with CBC. In 
Fig. 5b details of the spin-correlations are shown, with 
the most probable hole configuration projected out of the 
ORBA ground-state. At A=0.0 a very good convergence 
in the ORBA method can be achieved: starting with ini- 
tial states where the holes are uniformly distributed or 
phase separated, a fast convergence leads to the stripe 
results of Fig. 5b. 
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FIG. 5. (a) Rung hole density (n(r)) vs rung index r cal- 
culated using ORBA on a 4x6 cluster with 4 holes, J=0.2, 
t'=0.0, t"=0.25, and A=0.0, using CBC with k y =0. (b) 
Spin-spin correlations for 4 mobile holes at their most proba- 
ble configuration in the 4 holes ORBA ground state of (a). 

Summarizing, it has been shown that the striped states 
recently identified in models for the cuprates survive the 
introduction of an Ising anisotropy, particularly in the 
extended t-J model. The important physics that stabi- 
lizes stripes appears to be the competition between spins 
that order antiferromagnetically and holes that need to 
modify the spin environment to improve their movement, 
leading to an interesting potential extension into 2D of 
the familiar ID spin-charge separation ideas. The ap- 
proach discussed here focuses on the small J/t limit, and 
it appears unrelated with others based on the large J/t 
phase separated region. The fine details of the spin back- 
ground do not seem important, and as a consequence 
doped stripes should be a general phenomenon in corre- 
lated electronic systems. 
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As it occurs in the Nagaoka state of the one-hole sector of 
the t-J model, which as the cluster grows is stabilized at 
smaller J/t, reaching the J/t=0.0 limit in the bulk EB . 
The results thus far are not in contradiction with our 
recent interpretation of stripe formation as related to the 
one- hole properties jl^] , since this simple argument holds 
as long as the overall tendency in the spin sector is toward 
an AF state. Once ferromagnetism competes for one hole, 
the state that induces stripes upon further doping is now 
an excited state. 

Note that in average all AF correlations (across holes 
and between NN spins) uniformly decrease with A. The 
maximum AF correlation at zero doping for A=0.0 is 
-0.25 and for A=1.0 is x-0.35. 
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